import numpy as np
from matplotlib import pyplot as plt

plt.clf()
pKa = 4.76
pH=np.linspace(0,14,100) #calcul de 100 valeurs de pH de 0 à 14

pourcentageAH=1/(1+10**(pH-pKa))*100
pourcentageA=100-pourcentageAH

plt.title("Diagramme de distribution de l'acide éthanoique et de l'ion éthanoate")
plt.xlabel("pH")
plt.ylabel("Pourcentage des especes chimiques (%)")
plt.grid(True)
plt.plot(pH,pourcentageAH, c='red', label="Pourcentage en CH3CO2H")
plt.plot(pH,pourcentageA, c='blue', label="Pourcentage en CH3CO2-")
plt.legend()
plt.show()